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The review of modern study of algebraic, geometric and differential properties of 
quaternionic (Q) numbers with their applications. Traditional and "tensor" formulation of 
Q-units with their possible representations are discussed and groups of Q-units transformations 
leaving Q-multiplication rule form- invariant are determined. A series of mathematical and 
physical applications is offered, among them use of Q-triads as a moveable frame, analysis 
of Q-spaces families, Q-formulation of Newtonian mechanics in arbitrary rotating frames, and 
realization of a Q-Relativity model comprising all effects of Special Relativity and admitting 
description of kinematics of non-inertial motion. A list of "Quaternionic Coincidences" is 
presented revealing surprising interconnection between basic relations of some physical theories 
and Q-numbers mathematics. 

Introduction 

The discovery of quaternionic (Q) numbers dated by 1843 is usually attributed 
to Hamilton [TJ but in the previous century Euler and Gauss made a contribution 
to mathematics of Q-type objects; moreover Rodriguez offered multiplication rule for 
elements of similar algebra [3J HJ E]. Active opposition of Gibbs and Heaviside to 
Hamilton's disciples gave a start to the modern vector algebra, and later to vector 
analysis, and quaternions practically ceased to be a tool of mathematical physics, despite 
of exclusive nature of their algebra confirmed by Frobenius theorem. At the beginning 
of 20 century last bastion of Q-numbers amateurs, "Association for the Promotion of the 
Study of Quaternions" , was ruined. The only reminiscence of once famous hypercomplex 
numbers was the set of Pauli matrices. Later on quaternions appeared incidentally as a 
mathematical mean for alternative description of already known physical models (UJ [7j 
or due to surprising simplicity and beauty they were used to solve rigid body cinematic 
problems [B| . An interest to quaternionic numbers essentially increased in last two decades 
when a new generation of theoreticians started feeling in quaternions deep potential yet 
undiscovered (e. g. j§] - (TT|). 

This work is an attempt to give more systematic overview of contemporary state of 
Q-number mathematics, its applications to physical theories and possible perspectives in 
this area. In the context some quite specific even surprising physical models, but worth 
to pay attention to, are shortly discussed. 

The review arranged as follows. In section 1 general relations of the quaternionic 
algebra are briefly described in the traditional hamiltonian formulation as well as in 
tensor-like format. Section 2 is devoted to description of structure of three "imaginary" 
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quaternionic units. In section 3 the elements of differential Q-geometry are given 
with examples of their mathematical application. Section 4 comprises Q-formulation of 
Newtonian mechanics in the rotating frames of reference. Quaternionic Relativity Theory 
with a number of cinematic relativistic effects is found in section 5. Section 6 contains 
the list of " Great Quaternion Coincidences" and final discussion. 

1. Algebra of quaternions 

Traditional approach 

According to Hamilton, a quaternion is a mathematical object of the form 

Q = a + bi + cj + dk, 

where a,b,c,d are real numbers, a is a coefficient at real unit "1", and i,j,k - three 
imaginary quaternion units. The multiplication rule for these units given by Hamilton 
and often used in literature is 

li = il = i, lj=jl=j, lk = kl = k, 
i 2 =j 2 = k 2 = -l, 
ij = -ji = k, jk = -kj = i, ki = -ik = j 

These very cumbersome equations mean, that Q-multiplication loses a commutativity. 

Q1Q2 ^ Q2Q1, 

so that a notion of the right and the left multiplication appears, but it remains associative. 

(QiQ2)Q3 = Qi(Q 2 Q 3 ). 

Two rather different algebraic parts are separated naturally in a quaternion, these once 
could be denoted as scalar: 

seal Q = a, 

and vector 

vect Q = bi + cj + dk. 

Addition (subtraction) of quaternions is performed by components, scalar and vector parts 
are added (subtracted) separately. With respect to addition the Q-algebra is commutative 
and associative. 

Further step is quaternion conjugation introduced similarly to that of the complex 
numbers 

Q = seal Q — vect Q = a — bi — cj — dk, 
modulus of a Q-number is defined as 

\Q\ = ^/qQ = Va 2 + 6 2 + c 2 + d 2 . 

This permit to formulate a quaternionic division being as multiplication "right" and "left" 

n _ Q1Q2 n _ Q2Q1 
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Definition of Q-modulus enhances the famous four squares identity 

IQiQ 2 | 2 = |Qi| 2 |Q 2 | 2 . 

Due to the properties mentioned above the Q-numbers form the algebra, which belongs to 
the elite group of four the so-called exclusive - "very good" - algebras: of real, complex, 
quaternionic numbers and the octonions (Frobenious and Horwits theorems of 1878-1898 

C3). 

Special attention should be paid to Q-units representations. In terms of Hamilton 
real unit is simply 1 while three imaginary units similarly to complex numbers algebra 
are denoted as i, j, k. Later a simple 2x2 matrices representation of these units was 
revealed 

.( 1 \ 

J = - 

This representation of course is not unique. Here is a simple example. If in the above 
expressions imaginary unit i of complex numbers is represented as 2 x 2 with real elements 







then three vector Q-units turn out to be represented by real 4x4 matrices. The procedure 
of the matrix rank duplication can obviously be continued further. 

"Tensor" form and representations 

If each Q-unit is endowed with its proper number (as components of a tensor) 
(ij,k) -> (qi,q 2 ,q 3 ) = q, k,j,k,l,m,n,... = 1,2,3, 
then quaternionic multiplication rule acquires compact form 

where Sk n and Eknj - respectively, 3-dimension (3D) symbols Kronecker and Levi-Chivita. 

It is easy to show that a number of the Q-units representations even only by 2 x 2 
matrices is infinite. Indeed for any 2x2 matrices with properties 

A = [ a b I . D = I " ' I . TrA = TrB = 0. 




c —a 

the first two Q-units can be constructed as follows 

A B 

til = T, r » ^2 



VdetA VdetB 
while the third one is 

AB 

q 3 = qiq2 = —j^^=^^= provided that Tr(AB) = 0. 
V det AdetB 

The scalar unit is always invariant: 

1 = 




4 



Transformations of Q-units and invariancy of the multiplication rule 
a. Spinor-type transformations 

If U is an operator changing at once all the units, and there is an inverse operator 
U^ 1 : UU^ 1 = E, then transformations 

q fe , = UqtU- 1 and 1' = U1U' 1 = El = 1 

retain the multiplication rule 

form-invariant 

qfe'qn' = UqJJ^Uc^U = USknU' 1 + E kn jU O^U' 1 = 8 kn + SfenjOj'- 

Such operator can be represented for example by 2 x 2 matrix 



c d 



or unimodular quaternion, 
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a + d I f a + d 



u = ^ + \ i- -s- q. 



where 



a-d 
C 2 

In general this transformation contains 3 independent complex parameter (or 6 real ones), 
then U G SL(2,C). In special case of only three real parameters, then U G 577(2). 

6. Vector type transformations 

Vector Q-units can be transformed by 3 x 3 matrix O k > n 

The requirement of Q-multiplication form-invariance forces the transformation matrix to 
be orthogonal and unimodular 

Ok'nOj' n = S kn =>- 0~ k \ = O k ' n , det = 1. 

This transformation in general has 6 independent real parameters, then O G SO (3, C). In 
the special case of three parameters O G 5*0(3,7?). Below a variant of representation of 
the transformation matrix O is given with x, y, z being arbitrary real or complex functions 



O 



/ n 2 2 ^ — 2/ 2 — z 2 +yz\/l— x 2 — z 2 xj/— z\/l— x 2 — z 2 ^/ 1— j/ 2 — z 2 N 

v 1— z^ 1— z- 3 

Vl— x 2 — z 2 ^/ 1— j/ 2 — z 2 — xyz — yVl— x 2 — z 2 — 12-^/ 1— y' 2 — z 2 

x l^I 2 l^z 2 



z 



y a/ 1 - y 2 - z 2 J 
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This matrix can be represented as a product of three irreducible multipliers 
/ 

= 



l—x 2 —z 2 
1-z 2 



VT- 



V 



o 



1—X 2 —Z 
1-2 2 





o 

1/ 



-z 



\ 



-z 2 
1 
z v / T 3 ^ / 







1 1 1 









1— y 2 — z 2 
1-z' 2 



VT- 







l—y 2 —z 2 
1-z 2 



after substitutions z = sinB, rr = — sin A cos B, y = — sin T cos B, where A, £?,r - are 
complex "angles", it takes the form 



/ cos A sin A \ 



O = 



\ 



sin A cos A 
1/ 



^ cos BO— sin B ^ 

1 
\ sinB cosB / 



/l 







\ 



cos T sin r 
\0 — sin r cosr ) 



o$o*ol 



If the angles are real: A = a, B = j3, T = 7, then this transformation is an ordinary vector 
rotation consisting of three simple rotations around numbered orthogonal axes: O =^ 
R,R = R^R^RJ- Correlation between related "spinor" and "vector" transformations is 
easily determined: 



O k , n = --TriUcikU' 1 ^), U 



1 — Ok'n<jk(ln 

nun' 



Q-geometry in three dimensional space 

Hamilton was the first to note that triad of Q-units behaves as three strictly tied unit 
vectors (with length i) initiating Cartesian coordinate system, somewhat exotic because 
of its " imaginarity" . Due to the fact the Q-triad in 3D-space (qi,q2,qs) will be called 
'quaternionic basis' (Q-basis). Now Q-units transformations have apparent geometrical 
sense of various rotations of the Q-basis. An example: a simple rotation by real angle a 
around axis # 3 

q' = #3 q- 

Notion of Q-basis helps to introduce 3D quaternionic vectors (Q- vectors), defined as 



a = afcqfc, 



here all its components are real. The most important property of Q-vector - is its 
invariancy with respect to vector transformations from the group SO(3,R) 

a' = afc'qfc' = ak'Rk'jQj = a-j<ij = a. 

The projection of Q-vector onto arbitrary coordinate axis (represented by any different 
Q-unit) can be found in two ways. First, if at least one set of projections of Q-vector 
and rotation matrices R n k> are known then projections of this vector on rotated axis are 
immediately found 

O-k' — O'nRnk'- 

The second approach is related to existence of internal structure of the Q-units; a brief 
analysis of it is given in the next section. 
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2. Structure of quaternionic "imaginary" units 

Eigenf unctions of Q-units [13] 

Each vector Q-unit can be thought of as operator, so eigenf unctions and eigenvalues 
problem can be formulated for it 

qip = Xip, y?q = fiip. 

The solution of this problem are the eigenvalues ("imaginary length" of Q-unit with 
division by parity) 

A = \i — ±z, 

and two sets of eigenfunctions(one for each parity), possible given by columns ^ ± and 
rows (p , being the functions of components q. 

Here is an example explicit form of eigenfunction: for the Q-unit represented by 
matrix 

i I a b 



q= "T 



c —a 



where T = a 2 + be ^ 0, b ^ 0, c ^ 0, its eigenf unctions are defined as 

<t> k = x(\ ±^), ^ = y[\_ 

\ +T±a 

where x, y are arbitrary complex factors. 

The freedom of components, arising in the calculations is reduced by convenient 
normalization condition 

^ ± = 1, 

while the eigenfunctions orthogonality (by parity) is an inherited property 

= 0. 

One can construct tensor products of eigenfunctions and obtain 2x2 matrices 

G ± = VrV, 

possessing a properties reciprocal with respect to the ones of vector q: 

detC = 0, TrC = l, 

whereas 

det q = 1, Tr q = 0. 

Matrix C is idempotent 

C n = C, 

and can be expressed throw their own unit Q-vector 



C 



± 1 



7 



When inversed the latter expression gives information about internal structure of Q-unit 

q = ±i(2C ± - 1) = ±i(2iP ± if ± - 1), 

which turns out to consist of a combination of its eigenfunctions and scalar units. 

Since each Q-unit has its own eigenfunctions the Q-triad as a whole possesses 
unique set of eigenfunctions Wf^^fk)}- There is an interesting algebraic observation 
concerning this set. Three Q-units are interrelated by obviously nonlinear combination - 
multiplication e. g. 

qs = qiq2, 

but it is easy to show that corresponding eigenfunctions depend on each other linearly: 

Q-eigenf unctions help to represent a spinor-type transformation of Q-units retaining Q- 
multiplication invariant in the familiar form 

so that the eigenfunctions can be regarded as a set of specific spinor functions, allowing in 
subject in general to SL(2C) transformations. Yet another mathematical observation 
should be noted: from pairs of eigenfunctions, belonging to different Q-units of one 
triad and having one parity, one can construct 24 scalar invariants SL(2C) group; these 
invariants are real or complex numbers, e. g.: 



0T2 = = 

Quaternionic eigenfunctions as projectors 
Eigenfunctions act on their own Q-basis as following 

v^qi^fi) = ±», v^fi^Vfi) = 0, y^qs^fi) = o, 

or in general 

vfkfin^fk) = ±^fcn (no summation by k). 

It looks like that eigenfunctions select a projection of the unit Q-vector, generating 
them. This idea is confirmed by an example of an action of eigenfunctions of one Q-basis 
onto the vectors of the rotated Q-basis 

(pf k) q n ^f k) = yf k) R n im<hni)f k) = ±iRn'k = ±i cos Z(q„', q fc ) (no summation by k), 

the result of the action is 'nearly' projection of Q-basis q' on q. It is convenient to denote 
precise projection as 

(Qn')fc = ^^fk^n'^fk) = cos Z(q„', c^) (no summation by k). 

It is now easy to formulate rule of calculation of projection of a Q-vector a onto 
arbitrary direction, defined by vector q^ (e. g. with help of eigenfunctions of positive 
parity) 

(a)t = -^fc'V^qfc'^) = dk'Rk'j = dj (no summation by j). 

Thus quaternionic eigenfunctions with their own interesting properties are more 
fundamental mathematical objects then Q-units and too can serve as useful tool for 
practical purposes such as computing projections of Q-vectors. 
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4. Differential Q-geometry 

Quaternionic connection 

If vectors of Q-basis are smooth functions of parameters q^^) (index £ enumerates 
parameters), then 

dqfc($) = ujj. k j<ljd<$>z, 

where an object uj^ kj is called quaternionic connection. Q-connection is antisymmetric 
in vector indices 

uj^ kj + ^ jk — 0, 
and has the following number of independent components 

N = Gp(p-l)/2, 

where G is an number of parameters and p = 3 - is a number of space dimensions. If G 
= 6 [a case of group SO(3,C)], then N = 18; if G = 3 [a case of group SO(3,R)], then N 
= 9. Q-connection can be calculated at least in three ways: 

using vectors of Q-basis ue kn = \ 

\<y<±>£ 

using matrices U from the group SL(2C) (general case) and special representation of 
constant Q- units = — io~k, where Ok - Pauli matrices 

and, finally, using matrices O from 50(3,(7) (in a general case) 

90 tin 

All the formulae of course provide same result. 

From the point of view of vector transformations a Q-connection is not a tensor. If 
qfc = Okp>q. p >, then transformed components of connection are expressed throw original 
ones with addition of inhomogeneous term 

dOk p i 

U £, kj — O kp iOj n iU^ p ' n i + Ojpi — . 

In 3D space Q-connectivity has clear geometrical and physical treatment as moveable 
Q-basis with behavior of Cartan 3-frame. Parameters of its ordinary rotations can depend 
on spatial coordinates = then d n qk = QnkjQj, then components of slightly 

modified Q-connection 

Qnkj = ^ kjdn^i 

have a sense of Ricci rotation coefficients. Parameters can also depend on the length of 
line of motion of the Q-basis or on the observer's time. Then <3>£ = $^(t), d t c\k — &>kj<ij, 
and components of Q-connection 



&kj = kjdt^t 



9 



became generalized angular velocities of rotations of the frame. 
The typical examples of Q-frames and Q-connection are 

a) Frene frame. For the smooth curve x^(s) defined in constant basis the Frene frame 
is represented by the triad q^, obeying the equations 

^qi = Ri(s)q 2 , ^q2 = -Ri(a)qx + i?//(s)q 3 , ^q 3 = -R H (s)q 2 , 
where the first and the second curvatures are 



Rj — VI12, Rn — ^23. 

b) Twisted straight line. For a given straight line x± = u, £?> = £3 = 0, one can 
construct a Q-basis associated with it so that one vector is tangent to the line. In this 
case Q-connection is not zero and represented the only component describing torsion (or 
rather twist) of the line about itself. 

qi = ~ { I \ " I • <b = -i I )■ ^23 = 





du 

here 7(11) is the angle, which is an arbitrary but smooth function of the line length. 

Quaternionic spaces 

Tangent Q-space [15] . It is known that on every N-dimensional different iable 
manifold Un with coordinates {y A } one can construct a tangent space T N with coordinates 
{X^} so that dX^ = g^ dy B , where - Lame coefficients. By an extra rotation 
one can construct a tangent Q-space T(U, q), with coordinates {xk}, k = 1,2,3, which 
associated with Q-frame vectors. 

dx k = h k ( A) dX {A) = hk^g^dy 13 , 

where hyA) are in general non-square matrices normalized by projectors of the basic space 
onto 3D one or vice versa. 

Proper quaternionic space itself U3 is defined as 3D-space, locally identical to 
own tangent space T(U3,q). The Q-space has the following basic features. Its 
Q-metric represented by vector part of the Q-multiplication rule qj-q^ = — Sj k + £jkn<ln 
is nonsymmetric, its antisymmetric part is Q-operator (matrix), so that every point 
U3 has internal quaternionic structure. Q-connection U3 can be: (i) proper (metric) 
Q n kj = u>£kjdn&£, for variable Q-basis it is always non zero, and (ii) affine (non-metric), 
independent from Q-basis. Q-torsion does not vanish in both cases, whereas Q-curvature 
rknab = d a Q bkn - d b Q akn + Vt ajn Vt hjk - Q bjk Q ajn for the metric Q-connection identically 
equals zero, but can be present in the space of affine Q-connection. 

Once Q-space is introduced, there appears a new field of investigation of differential 
manifolds and spaces. Thus in the preliminary classification of Q-spaces based on presence 
and nature of curvature, torsion and non-metricity at least 10 different families can be 
distinguish 15J. In addition Q-spaces can be a nontrivial background for classical and 
quantum theories and problems. 
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4. Newton mechanics in Q-basis 

Dynamics equations in rotating frame ^B] 

The Q-basis endowed with clock becomes a classical (non-relativistic) reference 
system. For an inertial observer the dynamic equations of classical mechanics can be 
written in constant Q-basis 

m ^l x k^k = F k^k- 

SO (3, i?)-invariance of two Q- vectors, the radius- vector r = x k q k and force F = F k q k 
allow to represent these equations in Q-vector form 



d 2 , 

m -^( x k c lk) = F k q k , or ^ra- 



in explicit form these equations possess enough complicated structure 

d 2 d d 

1Ti\ ^£2^ n ^^ff/^k^kn ~\~ •^k'^^kn ^k^kj^jn) Fn 

which nevertheless can be simplified and interpreted from physical points of view. Due 
to antisymmetry of the connection (generalized angular velocity) 



f2j — fifcn c^knji ^kn ^j^-knji 



the dynamic equations can be rewritten in vector components 

m(a n + 2v k Q j e knj + x k -^e knj + x k n^ m e ]kv e mpn ) = F n 

or by conventional vector notation 

m{a + 2Q x v +fl x r + Q x (D, x f)) = F. 

Among left hand side terms one easily recognizes 4 classical accelerations: linear, Coriolis, 
angular and centripetal. However this traditional interpretation is good only for simple 
rotation; in the case of combination of many Q-frame rotations number of components 
of generalized accelerations highly increases, and the equations become much more 
complicated. However it is worth noting that derivation of these equations for the most 
complicated rotations with the help of Q-basis and Q-connection is extremely simple. 

Samples of Q-formulation of problems of classical mechanics 

'Chasing' Q-basis - is a frame with one of its vectors, say qi is always directed to 
observed particle. Dynamic equations for this case are written in explicit form in following 
manner 

f - r(tt 2 2 + Ql) = Fx/m, 
2rQ 3 + rQ 3 + rQ 2 ^i = F 2 /m, 
2rVL 2 + r(l 2 + rfiifi 3 = -F 3 /m. 
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Components of Q-connection are denned as functions of angles of two rotations, the first 
(an angle a) - around vector q 3 , the second (an angle 0) - around q 2 

VL\ = a sin f3, Q 2 = —$, ^3 = o- cos f3. 

The chasing Q-basis approach is convenient to solve a number of mechanical problems 
related to rotations, some times very complicated, of observed objects and systems of 
reference. Here is an illustration. 

Rotating oscillator. One seeks for motion law r(t) of a harmonic oscillator (mass m, 
spring elasticity k) which has a freedom of motion along rigid smooth rod rotating in the 
plane around one of its ends (here one end of the spring is fixed) with angular velocity 
cu] the equilibrium point is located at the distance 1 from the rotation center, there is no 
gravity. Radial and tangent dynamic equations in the chasing Q-basis (F is unknown rod 
reaction force) 

f — ruj = (r — 1), 2fuj = — F, 

m m 

admit the following family of solutions: 

(i) r(t) = r + v t + at 2 

mass moves away from the center of rotation with quadratic (or linear) law, 

(ii) r(t) = const + Ae iwt + Be~ iw \ w = yjk/m-uj 2 

here are three different situations depending on a relation of the quantities under the 
square root: 

- r = const, 

- harmonic oscillators, 

- exponential motion away from the center of rotation. 

It is interesting that the variants of rotating classical oscillator behavior with I = 
are precisely similar to behavior of four known cosmological models of Einstein-DeSitter- 
Friedman considered in the General Relativity. 



5. Construction of Quaternionic Relativity 



Hyperbolic rotations and biquaternions ^7] 

It was noted above, that SO(3,C)-transformations of Q- units admit pure imaginary 
parameters. In this case rotations become hyperbolical (H - from hyperbolic); e. g. 
simple H-rotation q' = H$ q is performed by matrix of the form 

' cosh ip —i sin ip ^ 
Hf = i sin if} cosh if> , 
V 1/ 

and 2x2 -matrices of Q-units representation are no longer hermitian: 







-i> 



12 



This is the time to recall the notion of so called biquaternionic vectors (BQ). BQ- vector 
is defined as Q- vector with complex components u = (a/. + ibk)q.k- Obviously for vectors 
of this type the norm (or modulus) in general can not be defined; but among all BQ- 
vectors there is a subset of "good" elements with well definable norm by u 2 = b 2 — a 2 . 
These vectors appear to be form-invariant with respect to transformations of subgroup 
SO(2, 1) C SO(3,C), and in particular, with respect to simple H-rotations q' = Hqu = 
Uk<ik = Uk'(\k', but only when reciprocally imaginary components dkbk = are orthogonal 
to each other. 

Quaternionic Relativity 

The made above observation allows to suggest a space-time BQ- vector "interval" 

dz = (dx k + idt k )qk, 

with specific properties: 

(i) Temporal interval is defined by imaginary vector, 

(ii) space-time of the model appears to have six- dimensional (6D), 

(iii) vector of the displacement of the particle and vector of corresponding time 
change must always be normal to each other dxkdtk = 0. 

In this case BQ- vector-interval is invariant under group £0(2, 1) C SO(3, C), as well 
as of course its square (which differs from the square of norm only by sign) dz 2 = dt 2 — dr 2 , 
the latter has precisely the same form as a space-time interval of Special Relativity of 
Einstein. This 6D-model was initially named the Quaternionic Relativity. Temporal and 
spatial variables symmetrically enter the expression of BQ- vector-interval, and the Q-triad 
related to it describes relativistic system of reference X = (qi,q2, Transition from 
one reference system to another is performed with the help of 'rotational equations' of 
the type £' = OS with matrix O from the group SO(2, 1) is a product of matrices of real 
and hyperbolic rotations. So the theory could also be named (may be more correctly) 
'Rotational Relativity'. The meaning of a simple H-rotation is immediately revealed from 
the first line of equation £' = HfT, in the explicit form 

iqx' — icosh?/>(qi + tanh^q 2 ). 

If like in Special Relativity cosh^ = dt/dt', then 

idt'qr = idt(qi + Vq 2 ), 

which describes motion of reference system £' relative to E with velocity V along direction 
q 2 . It is easy to show that SO(2, l)-rotations of Q-reference system enhance Lorenz 
coordinate transformations and therefore all cinematic effects of Special Relativity. 

It should be noted here that parameters of real and hyperbolic rotations can be 
variable for instance dependent on observer's time. This hints to expect of the discussed 
theory a possibility to describe non-inertial motions. Analysis of the rotational equations 
confirms the expectation. Well-known relativistic model of reference system constantly 
accelerated with respect to the inertial one (hyperbolic motion), frequently found in 
literature and normally regarded with use of assumption beyond frames of Special 
Relativity, in quaternionic theory is solved naturally and fast not only from the inertial 
observer viewpoint, but from position of accelerated frame too [TH] . 
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The kinematic problem of other non-inertial motion - relativistic circular motion 
can be completely and precisely resolved by means of the rotation equation 
£' = H^^Ri^Ti, where £' is reference system rotating along the circle around the 
immobile frame E. This problem also can be solved both from the point of view of 
inertial observer, in this case the result has the form 

dt' cosh ip(f), a ^ = ^/ dt ' tanh -^P( t ')^ 

1 dip , . / da(t)\ 2 

ftanl^J 7~o ', 77 1 Q"norm\J> I -K I 7 I ) 

cosh 2 ^ dt \ dt J 

and from the point of view of the observer in the reference system arbitrary moving along 
circular orbit. 

The solution of the problem of "classical" Thomas precession in the framework of 
Special Relativity also needs additional assumptions, while in the quaternionic theory has 
a single line form - the first row of the matrix of rotation equation S" = R^ a ^ -R" E, 
in this case of course correct value of precession frequency is obtained 

ujt = (1 — cosh-?/;) rs — -luV 2 . 

Moreover, the quaternionic theory of relativity appears to be able to describe Thomas 
precession for the vectors moving along trajectories of general type. The basic rotational 
equation in this case naturally generalized: E" = R~ 6 ^' //w) R 6 ^>Y1, here 9(t) - an angle of 
instant rotation. Requirement that an axis of hyperbolic rotation be normal to the plane 
formed by the radius- vector of observed frame and its velocity vector, is also significant. 
In this case formula of variable frequency of general Thomas precession has the form 

n, = !(»-*). 

An example of such Thomas precession is an apparent displacement of mercurial 
perihelion, for which calculations give a value Ae = 2, 7"/100 years. 

Universal character of motion of the bodies (including non-inertial motions) in the 
Quaternionic Relativity suggests seeking for new cinematic relativistic effects. One is 
found in Solar System planets' satellites motion. Relative velocity of the Earth and other 
planets changes with time and sometimes achieves significant value comparable somehow 
to value of the fundamental velocity. This can lead to discrepancy between calculated 
and observed from the Earth cinematic magnitudes characterizing cyclic processes on 
this planet or near it. In particular there must be a deviation of the planetary satellite 
position. Such an angular difference is surprisingly found to be linearly dependent upon 
the time of observation 

luVeVp 
« —t, 

c 2 

here u is an angular velocity of satellite motion around the planet, V - are linear velocities 
of the Earth and the planet around the sun. The magnitude of the effect is the following 
for the closest to the Jupiter and "the fasters" Jupiter satellite Aip = 12' for 100 terrestrial 
years; for the Mars satellite (Phobos) A(p = 20' for 100 terrestrial years JH] • Both values 
are big enough for the effect to be noticed in prolonged and precise observations. 
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One can say that space-time model and kinematics of the Quaternionic Relativity are 
nowadays studied in enough details and can be used as an effective mathematical tool for 
calculation of many relativistic effects. But respective relativistic dynamic has not been 
yet formulated, there are no quaternionic field theory; Q-gravitation, electromagnetism, 
weak and strong interactions are still remote projects. However, there is a hope that it 
is only beginning of a long way, and the theory will "mature". This hope is supported 
by observation of number of remarkable "Quaternionic Coincidences" forming a discrete 
mosaic of physical and mathematical facts; probably one day it will turn into a logically 
consistent picture providing new instruments and extending our insight of physical laws. 

6. Remarkable "quaternionic coincidences" 

There are, at least, five such coincidences (all of them given below), noted by different 
authors in various time. 

1. The Maxwell equations as an conditions of the analyticity of functions of 
quaternionic variable. 

In 1937 year Fueter [2II| noted, that Cauchy-Riemann df/dz* = equations defining 
the differentiability of complex variable function and modeling physically a flat motion of 
liquid without sources and whirls, have the following quaternionic analogue 



Surprising fact is that the equations of classic Maxwell electrodynamics in vacuum prove 
to be corresponding physical model 



2. Classical mechanics in the rotating reference systems. 

The compact form of Newton equations in quaternion frame is described above in 
section 4. Finally it should be stressed that the form of dynamics equations naturally 
arising and externally primitive 



hides all possible combinations of rotations of reference systems or observed bodies. Using 
differential quaternionic objects helps to easily obtain explicit form of the equations whose 
elements have obvious physical meaning. 

3. The quaternionic theory of relativity. 

1:1 isomorphism of the Lorenz group of Special Relativity and the group of invariance 
of quaternionic multiplication 5*0(3,(7) leads to non-standard theory of relativity with 
symmetric six-dimensional space-time. This theory significantly differs from Einstein 
Special Relativity in origin, model, possibilities and mathematical tools, but predicts 
absolutely similar cinematic effects. Invariance of specific biquaternionic vector "interval" 
dz = (dxkn + i dtk)c[k under subgroup 5*0(2, 1) with in general variable parameters admits 
calculation of relativistic effects for non-inertial motion of reference systems. 




divE = 0, divB = 0, rotE - — = 0, rotB + — = 0. 



dt dt 



mr = F 
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4. Pauli equations [2T| . 

Consider the quantum particle with electric charge e, mass m, and generalized 
momentum 

€ 

Pk = -ih- A k 

OXk c 

in the simplest quaternionic space (all the parameters are constant, connection, non- 
metricity, torsion and curvature equal to zero). Hamiltonian of such particle in Q-metrics 

H = — PkPm Qfe^m 

2m 

is the exact copy of Hamilton function of Pauli equation 

H=±(p- e -A) 2 -^B.a, 

2m V c / 2mc 

and the spin term "automatically" acquires a coefficient equal to Bohr magneton. 

5. Young-Mills field strength. 

If one constructs a "potential" vector in an arbitrary quaternionic space from 
Q-connection components VL amn (indices a, b, c enumerate coordinates of basic Q-space, 
indices j, k, m, n enumerate vectors of tangent triad) 

A-ka = '^■kmrS^amni 

and similarly construct a "field strength" vector 

_ 1 

Fkab = -e kmn r mnab , 
from quaternionic curvature components 

Tknab ^a^&fcn db^akn ^ajrS^bjk ^bjk^ajn 

then these two geometrical objects are interconnected in the similar manner as the field 
strength and potential of the Young-Mills field 

Fkab = 9 A ka 9 a j4^b £kmnA- ma A n b. 

(formula) It should be stressed that for the Q-spaces with metric (not affme) connection 
curvature (field strength) identically vanish. 



Discussion 

Quaternionic numbers of course are first of all mathematical objects, so the problem 
of development of their algebra, analysis and geometry is self-consistent. But history 
of modern science states that once the geometry, in particular differential geometry, is 
discussed the presence of physics is unavoidable. There is a known point of view that 
Einstein who suggested General Relativity was a pioneer in geometrization of physics. 
But it is also known that quite earlier Maxwell formulated his electrodynamics in terms of 
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quaternions convenient for description of 'etheric tensions' which were thought to represent 
field strength vectors. But since that the geometrical language has not been utilized for 
many decades. 

The aspects of quaternionic mathematics given in this review once again draw 
attention to 'genetic relations' between physics and geometry: from description of frames 
rotations to quaternionic field structure phenomena in Pauli equations and Young-Mills 
theory. 

Wide variety of possibilities provided by Q-approach and derived within it 
non-traditional physical models, like six-dimensional space-time or mentioned above 
coincidences may lead to opinion that quaternions are still a mathematical play, something 
like 'lego' elements, from which one can build many exotic constructions. 

As a comment there are the following two observations. 

1. Producing non-standard physical models Q-method nonetheless allows to 
successfully solve physical problems thus being a useful tool for practical purposes. A 
typical example: inherited exponential character of representation of simple rotations 
helps to simply formulate summation of different rotations, including, of course, imaginary 
rotations, describing relativistic boosts. Recall that in classical mechanics summation of 
ordinary rotations is quite a task. 

2. All physical quaternionic theories are not heuristically invented, but appear 
naturally from fundamental mathematical lows, as though confirming Pythagorean idea 
on "world - number" dependence. Indeed, Q-algebra, the last associative algebra, 
describes well physical quantities, all of them up to our knowledge being associative 
with respect to multiplication: from observable kinematic and dynamic one, to tensors 
and spinors incorporated in the theories. All this gives a hope that further efforts in 
the research "quaternions - physical laws" relations will once grow into wide scientific 
programme. Yet another small, but persevering step in this direction has been recently 
made, when the author of this review succeeded to found an exact solution for relativistic 
oscillator problem in the framework Quaternionic Relativity. Details of the solution will 
be published elsewhere. 
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